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1.  INTRODUCTION 


1.1  BACKGROUND  AND  OBJECTIVES 

In  the  initial  phase  of  this  study,  a  range  and  bear¬ 
ing  (RB)- depen dent  model  was  developed  to  predict  the  time 
difference  (TD)  grid  for  the  St.  Marys  River  Loran-C  minichain 
(Ref.  1).  The  RB  model  was  calibrated  using  TD  data  collected 
over  an  18-day  period  in  September  and  October  1977*,  and 
therefore  exhibits  the  "average"  propagation  conditions  charac¬ 
terized  by  the  data  at  that  time.  Since  the  chain  is  expected 
to  provide  an  all-weather  navigation  capability,  the  sensitivity 
of  the  TD  grid  to  seasonal  variations  in  propagation  conditions 
is  of  much  interest.  In  Ref.  1,  a  comparison  of  the  data 
collected  by  the  U.S.  Coast  Guard  in  November  1977  to  that 
collected  in  September  showed  rather  large  (up  to  450  nsec)  TD 
variations  at  some  data  collection  sites.  These  observed  vari¬ 
ations  provided  additional  motivation  for  the  present  study, 
which  presents  a  theoretical  sensitivity  analysis  of  the  TD 
grid  based  on  expected  seasonal  variations  in  propagation  con¬ 
ditions. 


1.2  TECHNICAL  APPROACH 

Figure  1.2-1  presents  an  overview  of  the  approach 
employed  to  predict  the  variations  in  the  St.  Marys  River  TD 
grid  induced  by  seasonal  changes  in  propagation  conditions. 
The  propagation  medium  of  the  Loran-C  signal  is  first  charac¬ 
terized  by  the  critical  signal  propagation  parameters.  After 

•This  period  is  hereafter  referred  to  as  September  1977. 
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Figure  1.2-1 


Sensitivity  Analysis  Approach 
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establishing  a  functional  relationship  between  the  parameters 
and  TDs,  sensitivity  of  the  grid  to  changes  in  propagation 
parameters  is  determined.  The  sensitivity  equations  reveal 
the  specific  role  played  by  the  System  Area  Monitor  (SAM)  in 
controlling  grid  sensitivity  (confirming  heuristic  arguments 
in  Ref.  1). 

Next,  expected  seasonal  variations  in  the  propagation 
parameters  are  derived  from  historical  meteorological  data 
and  other  published  information.  Each  parameter  is  charac¬ 
terized  by  a  mean  seasonal  cycle  and  by  random  fluctuations 
about  this  mean  cycle  encountered  over  many  years.  By  employ¬ 
ing  this  information  in  the  sensitivity  equations,  mean  and 
random  seasonal  TD  variations  are  predicted  for  the  four  sea¬ 
sons.  A  comparison  of  predicted  TD  variations  with  recorded 
seasonal  data  shows  that  propagation  theory  carnot  be  used  to 
explain  all  the  observed  variations.  Preliminary  recommenda¬ 
tions  for  a  comprehensive  data  acquisition  program  are  pro¬ 
vided,  based  on  the  sensitivity  analysis. 

Position  error  statistics  at  St.  Marys;  River  waypoint 
are  computed  by  interpreting  seasonal  TD  variat ions  as  uncom¬ 
pensated  errors  in  the  calibrated  grid.  These  statistics  are 
compared  with  similar  statistics  presented  in  lief.  1  to  deter¬ 
mine  the  relative  contribution  of  grid  instabi'.  ity  to  position 
errors . 


1.3  REPORT  OVERVIEW 


The  theoretical  framework  employed  to  predict  TD  vari 
ations  —  including;  definitions  of  propagation  parameters,  the 
relation  between  parameters  and  TDs,  and  an  analysis  of  the 
effect  of  the  SAM  ---  is  included  in  Section  2.  In  Section  3. 
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seasonal  variations  in  propagation  parameters  are  inferred 
from  available  information,  and  TD  and  position  errors 
resulting  from  the  variations  are  determined.  Conclusions 
and  recommendations  of  the  study  are  summarized  in  Section  4. 
Detailed  equations  and  tables  are  presented  in  the  appendices. 


THE  ANALYTIC  SCIENCES  CORPORATION 


2.  PROPAGATION  MODELS  AND  SENSITIVITY  ANALYSIS 

2 . 1  INTRODUCTION 

The  range  and  bearing-dependent  model  calibrated 
in  Ref.  1  is  a  semi -empirical  characterization  of  the  St.  Marys 
River  Loran-C  time  difference  grid.  Although  the  functional 
structure  of  the  model  is  based  on  classical  propagation  theory 
(Ref.  2),  it  does  not  explicitly  include  the  propagation  para¬ 
meters.  The  model  coefficients  were  chosen  to  minimize  the 
residuals  between  predicted  and  measured  TDs ,  thereby  charac¬ 
terizing  the  grid  which  existed  during  the  September  1977  data 
collection  period.  Since  the  model  does  not  explicitly  include 
all  important  propagation  parameters,  it  cannot  be  used  directly 
in  a  seasonal  sensitivity  analysis. 

Two  alternate  approaches  to  seasonal  sensitivity 
analysis  have  been  considered: 

•  Calibrate  a  t ime-dependent  TD  prediction 
model  which  includes  the  important  pro¬ 
pagation  parameters,  using  seasonal  TD 
data 

•  Utilize  a  model  of  a  homogeneous  smooth 
earth  to  estimate  gross  changes  in  the 
TD  grid  due  to  uniform  changes  in  pro¬ 
pagation  parameters. 

The  former  approach  would  be  capable  of  accounting  for  spatial 
detail  of  the  ground  conductivity,  through  a  bearing  dependence 
as  in  the  RB  model  or  through  some  similar  artifice.  However, 
this  requires  a  more  extensive  data  base  than  is  presently 
available.  The  latter  approach,  although  not  addressing  changes 
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in  spatial  detail  explicitly,  has  been  chosen  here  5  or  its 
tractability .  The  model  is  not  calibrated  against  seasonal 
data  but  rather  stands  on  theory  alone . 


2.2  ACCURACY  OF  CALIBRATED  TD  GRID 

The  role  of  the  RB  model  in  predicting  the  St.  Marys 
River  Loran-C  TD  grid  is  reviewed  in  this  section.  The  true 
time  difference  (TD.^),  between  signals  from  secondary  station 
i  and  the  master  station  m,  is  assumed  to  be  measured  by  a 
perfect  receiver  in  a  noiseless  environment  (this  ideal  con¬ 
dition  is  referred  to  as  zero  measurement  noise).  If  propa¬ 
gation  delays  between  the  stations  and  the  user  are  denoted 

by  PD.  and  PD  and  the  secondary  station  emission  delay  by 
1  m 

ED^ ,  then 


TD.  *  PD.  -  PD_  +  ED.  (2.2-1) 

x  x  m  x 

All  variables  in  Eq .  2.2-1  may  vary  temporally  at  a  fixed  user 
location . 


In  Ref.  1,  the  true  time  difference  is  conceptually 
separated  into  three  components: 

•  A  portion  which  can  be  described  by 
the  RB  model  structure,  termed  the 
spatial-area  TD,  TD^ 

•  A  portion  due  to  unmodelled  spatial 
details  which  do  not  change  with  time, 
termed  the  local  warpage ,  6TDi 

•>  A  portion  due  to  temporal  variations 
in  propagation  conditions,  termed 
the  temporal  grid  instability,  ATD^ 

The  RB  model  predicts  tie  spatial-area  TPs  which  existed  during 
September  1977  data  collection  (RB  model -predicted  TD ,  denoted 

A 

by  TDi).  The  error  incurred  in  estimating  the  spatial-area  TD 
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(termed  residual  grid  prediction  error  or  6TTL)  depends  on 
the  quality  and  quantity  of  data  used  to  calibrate  the  RB 
model.  This  error  was  predicted  at  the  waypoints  by  covari¬ 
ance  analysis  in  Ref.  1.  The  total  TD  prediction  error 
(difference  between  the  predicted  and  measured  TD,  in  the 
absence  of  measurement  noise)  when  the  RB  model  is  employed 
to  predict  the  true  time  difference  is  then 

/v 

Total  TD  error  =  TDi  -  TD±  =  STIK  +  6TDi  +  ATD±  (2.2-2) 

The  primary  purpose  of  this  study  is  to  predict  sea- 
onal  statistics  for  temporal  grid  instability.  However,  during 
the  course  of  computing  total  position  errors,  local  warpage 
errors  are  hypothesized  based  on  available  information. 


2.3  PROPAGATION  PARAMETERS  AND  THEORETICAL  PIASE 
DELAY  MODEL 

Signal  propagation  theory  (Ref.  2)  and  published 
observations  (Refs.  4  and  5)  show  that  the  characteristics 
of  Loran-C  groundwave  propagation  can  be  described  by  the 
following  propagation  parameters: 

•  Surface  atmospheric  refractive  index,  n 

•  A  parameter  a,  which  is  related  to  the 
refractive  index  gradient 

•  Ground  conductivity,  a  . 

The  parameters  vary  both  spatially  and  temporally  in  the 
St.  Marys  River  Loran-C  coverage  area.  However,  the  classi¬ 
cal  propagation  model  used  herein  assumes  a  homogeneous 
smooth  earth,  an  idealization  of  the  true  propagation  medium. 
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Classical  theory  for  a  homogeneous  medium  expresses 
the  groundwave  signal  phase  delay*  as  the  sum  of  a  primary 
phase  delay  (T)  and  a  secondary  phase  delay,  also  known  as 
secondary  phase  factor,  SF+.  The  primary  phase  delay  is 
experienced  in  an  atmosphere  with  refractive  index  n,  in 
the  absence  of  the  earth,  while  the  secondary  phase  delay 
describes  the  additional  effect  of  the  earth.  The  theoreti¬ 
cal  time  difference  is  therefore  given  by 

TD±  =  (T.-Tm)  +  (SFi-SFm)  +  ED±  (2.3-1) 

In  Sections  2.3.1  to  2.3.3,  the  relationships  between  the  pro¬ 
pagation  parameters  and  the  components  of  phase  delay  are 
presented.  Then,  in  Section  2.4,  TD  sensitivity  to  changes 
in  the  parameters  is  derived  from  Eq.  2.3-1. 

2.3.1  Surface  Refractive  Index 

The  wave  velocity  in  an  atmosphere  of  refractive  in¬ 
dex  n  (dimensionless)  is  reduced  from  the  velocity  in  free 
space  (c  in  ft/ysec)  by  the  factor  1/n.  Hence,  the  primary 
phase  delay  (ysec)  along  a  geodesic  of  length  R  (ft)  is 

T  =  (2.3-2) 

Secondary  phase  delay  is  approximately  proportional  to  n1^ 

(Ref.  2).  Therefore,  changes  in  SF  are  given  by 

ASF  -  f  (x)  <2-3-3> 


♦Phase  delay  $  is  normally  expressed  in  units  of  time  through  the 
conversion  4>/2Trf,  where  f  is  the  100  Khz  Loran-C  carrier  fre¬ 
quency  . 

tHerein  SF  denotes  the  total  secondary  phase  factor  for  any 
land/water  path  and  not  the  SF  associated  with  an  equivalent 
sea  water  path. 
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Since  changes  in  primary  phase  delay  (from  Eq.  2.3-2)  are 
given  by 


AT  =  T  ~  (2.3-4) 

and  T>>SF/3  for  ranges  of  interest  in  the  St.  Marys  River  cover¬ 
age  area,  the  effect  of  changes  in  the  refractive  index  on  SF 
can  be  ignored. 

2.3.2  Parameter  a 

The  parameter  a  (dimensionless)  is  an  artifice  used  in 
signal  propagation  theory  to  account  for  wave  refraction  due  to 
the  refractive  index  gradient,  dn/dh  (where  h  is  altitude).  For 
an  earth  of  radius  a,  the  parameter  is  computed  by  evaluating 

“  =  1  +  I  3E  <2'3-5> 

at  the  earth's  surface  (h=0).  The  parameter  a  affects  second¬ 
ary  phase  delay  in  a  complex  manner  (see  Ref.  2  for  details). 
However,  over  ranges  of  interest  in  the  St.  Marys  River 
Loran-C  coverage  area  (<100  mi),  the  following  relation  pro¬ 
vides  a  reasonable  approximation  to  the  classical  model: 

SF(a )  =  SF(aQ)  +  0.0005  (a-aQ)  T  (ysec)  (2.3-6) 

where  aQ  is  a  nominal  value  of  a  (chosen  equal  to  0.75,  to 
lie  in  the  expected  range  for  the  St.  Marys  River  region)  and 
T  is  defined  by  Eq .  2.3-2. 

2.3.3  Ground  Conductivity 

Ground  conductivity  c (mho/m)  is  the  most  important 
parameter  in  the  theoretical  formulation  of  secondary  phase 
delay,  encompassing  all  non-atmospheric  effects.  The  theore- 
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tical  relationship  between  secondary  phase  delay  and  conduc¬ 
tivity,  detailed  in  Ref.  2,  is  too  cumbersome  to  apply  in 
sensitivity  analyses.  Therefore,  an  alternative  relationship 
has  been  derived  which  is  only  applicable  to  the  signal  path 
lengths  (1-100  mi)  and  conductivities  (0.001-0.01  mho/m,  see 
Section  3.2.3)  relevant  to  the  St.  Marys  River  region.  The 
following  polynomial  was  fit  to  "samples"  of  secondary  phase 
delay  computed  from  classical  theory: 

2  1  # 

SF  -  £  £  V  Tk  (2.3-7) 

k=-2  £=- 1 

where 

SF  =  secondary  phase  delay  (ysec) 

T  =  primary  phase  delay  (ysec) 

5  =  natural  logarithm  of  conductivity,  oc  (mho/m) 

A,  =  least-squares-derived  coefficients  listed 
in  Table  2.3-1. 


TABLE  2.3-1 

COEFFICIENTS  A,  „  OF  POLYNOMIAL  SF  MODEL 

kJl 


X 

-1 

0 

1 

-2 

-106.2 

-58.80 

-9.115 

-1 

33.12 

21.18 

2.850 

0 

-3.766 

-2.139 

-0.3279 

1 

-0.02380 

-0.01380 

-0.002257 

2 

1 . 565xl0-5 

9. 369xl0-6 

1 . 456xl0~6 
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Since  the  coefficients  AkJl  were  computed  based  on  the  nominal 
value. of  a  (aQ  =  0.75),  the  secondary  phase  delay  defined  by 
Eq.  2.3-7  is  the  nominal  delay,  SF(aQ),  appearing  in  Eq.  2.3-6 
Figure  2.3-1  illustrates  the  relationship  between  SF  and  T, 
computed  from  Eq .  2.3-7,  for  two  pertinent  values  of  conduc¬ 
tivity.  Polynomials  similar  to  Eq .  2.3-7  have  been  employed 
in  other  Loran-C  studies,  such  as  described  in  Ref.  6. 


M47N0 


Figure  2.3-1  Secondary  Phase  Delay  as  a  Function 

of  Primary  Phase  Delay  for  Two 
Conductivities 

2.4  SENSITIVITY  EQUATIONS  AND  THE  EFFECT  OF  SAM 

In  Section  2.3,  the  relationship  between  TDs  and  the 
propagation  parameters  of  a  homogeneous  medium  #as  established 
Before  Eq.  2.3-1  can  be  applied  to  determine  the  sensitivity 
of  TDs  to  temporal  variations  in  the  propagation  parameters, 
the  secondary  station  emission  delay  (see  Eq.  2.2-1)  must  be 
related  to  the  parameters. 
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Emission  delay,  ED^  is  the  time  between  transmission 
of  the  master  pulse  group  and  the  secondary  pulse  group.  The 
System  Area  Monitor  (SAM)  controls  emission  delay  through  cod¬ 
ing  delay  adjustments  transmitted  to  the  secondary  station. 

These  adjustments  are  chosen  so  that  the  TD  measured  at  SAM 
equals  a  prescribed  reference  (TD?)  with  a  tolerance  of  +15  nsec. 
If  SAM  were  capable  of  controlling  its  received  TDs  precisely , 
then  from  Eq.  2.3-1 

EDi  =  TD?  -  (T.'-T^)  -  (SF!-SF^)  (2.4-1) 

where  prime  (')  distinguishes  SAM  phase  delays  from  user  phase 
delays.  Note  that  the  emission  delay  varies  if  the  propagation 
parameters  of  the  medium  vary.  By  combining  Eqs .  2.3-1  and 
2.4-1,  the  effect  of  propagation  parameters  and  SAM  :>n  the  TDs 
at  a  user  location  is  shown  to  be 

TD.  -  f*  (T .  -T  )  +  (SF.-SFjl  -  f(Ti-T')  +  (SF'-SF')l  +  TD°  (2.4-2) 

l  L.  *  ®  l  m  j  l  l  m  l  m  j  1 

Note  that  a  user  positioned  at  SAM  measures  TD°. 

Equation  2.4-2  can  be  employed  directly  in  sensitivity 
analyses  to  predict  ATD^^ ,  the  temporal  grid  instability  (see 
Section  2.2)  arising  from  a  particular  temporal  change  in  propa¬ 
gation  parameters.  If  subscripts  1  and  2  denote  two  different 
parameter  sets,  then 

ATD.  =  TD.(n2,a2,£2)  -  TD^nj  ,0^  ,  q)  (2.4-3) 

If  random  changes  are  of  interest,  it  is  more  convenient  to  use 
the  following  linear  approximation  in  a  covariance  analysis: 

3TD .  3TD .  3TD. 

ATDi  -  IT  An  +  T 3T  Aa  +‘  —  (2.4-4) 
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The  required  partial  derivatives  of  TD  with  respect  to  the 
parameters  follow  from  Eq.  2.4-2,  using  the  relations 


=  0.0005  T 
da 

!r  -  k?_2  ***  1 <2-*-5 

Figure  2.4-1  shows  the  variation  of  the  derivative  of  SF  with 
respect  to  £  with  T,  for  two  values  of  conductivity  referred  to 
later  herein.  The  convenient  approximation 


9SF 

3£ 


-0.039  T0-47 


(2.4-6) 


differs  from  the  plotted  values  by  less  than  five  percent  for 
both  values  of  conductivity. 


figure  2.4-1  Sensitivity  of  Secondary  Phase 

Delay  to  fcn  (Conductivity) 
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Equations  2.4-2  to  2.4-6  can  be  used  to  relate  param¬ 
eter  variations  to  expected  TD  changes  and  are  referred  to  as 
the  sensitivity  equations.  The  following  observations  can  be 
made  from  the  sensitivity  equations: 

•  The  reference  TD  prescribed  for  SAM 
does  not  affect  sensitivity 

•  The  SAM  removes  all  sensitivity  of 
TDs  to  the  relative  drift  of  the 
master  and  secondary  clocks 

•  SAM  decreases  the  sensitivity  of 
TDs  to  propagation  parameter  varia¬ 
tions  at  some  user  locations,  while 
increasing  sensitivity  at  others.* 

The  impact  of  SAM  on  TD  sensitivity  to  variations  in  refractive 

index  and  the  parameter  a  is  determined  by  the  primary  phase 

delay  differences  at  the  user  (T.-T_)  and  at  the  SAM  (T!-T_L). 
- * -  '  i  m  l  m 

SAM  increases  the  magnitude  of  the  TD  sensitivity  to  n  and  a 
at  those  user  locations  where 

•  The  primary  phase  differences  at  the 
user  and  SAM  are  of  opposite  sign 

•  The  primary  phase  differences  are  the 
same  sign  but  | T ^ — Tm |  <  0. 5 j T^-T^ | . 

SAM  decreases  the  TD  sensitivity  at  all  other  locations  in  the 
coverage  area,  reducing  the  sensitivity  to  zero  alon~,r  the  SAM 
LOP.  In  other  words,  because  SAM  is  controlling  the  chain  to 
hold  the  LOPs  through  SAM  constant,  TD  variations  due  to  pro¬ 
pagation  parameter  changes  are  removed  by  SAM  and  a  user  on 
the  SAM  LOP  will  not  observe  the  effects  of  these  parameter 
variations  (within  the  tl5  nsec  control  tolerance).  Conse¬ 
quently,  the  associated  TD  sensitivity  to  parameter  changes 
is  zero.  However,  there  is  a  locus  of  points  in  the  coverage 
region,  for  each  TD,  along  which  the  magnitude  of  the  TD 


*This  conclusion  was  arrived  at  in  Ref.  1  by  assuming  a 
linear  dependence  of  phase  delay  on  a  propagation  param¬ 
eter. 
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sensitivity  is  the  same  with  SAM  as  without  SAM.  This  is  not 
to  say  that  the  TD  sensitivity  is  zero  and  that  a  user  will 
not  observe  TD  variations  along  the  locus  of  points.  In  fact, 
the  sense  of  the  TD  sensitivity  changes  if  SAM  is  active  or 
inactive  but  the  absolute  value  of  the  sensitivity  remains  con¬ 
stant.  In  Fi*g.  2.4-2,  the  St.  Marys  River  region  is  divided 
into  regions  (for  each  TD)  where  SAM  increases  or  decreases 
the  TD  sensitivity  to  variations  in  n  and  a.  From  these  charts 
it  is  evident  that  SAM  increases  grid  sensitivity  to  variations 
in  n  and  a  at  many  waypoints  in  the  southern  region  of  the 
river. 


The  sensitivity  of  TDs  to  temporal  variations  in  the 
conductivity  of  a  homogeneous  medium  is  nonlinear  in  T,  as 
shown  by  Eq.  2.4-6.  Therefore,  the  boundary  lines  in  Fig.  2.4-2 
are  not  precisely  applicable  to  conductivity.  In  Fig.  2.4-3, 
the  sensitivity  of  TDs  to  conductivity  variations  is  plotted 
versus  waypoint,  for  cases  with  SAM  and  without  SAM.  The  fol¬ 
lowing  observations  are  derived  from  these  sensitivity  curves: 

•  TDX  is  the  least  sensitive  TD.  Although 
SAM  increases  the  sensitivity  to  conduc¬ 
tivity  variations  at  most  waypoints, 
sensitivity  is  decreased  to  near  zero 
between  waypoints  I  and  M,  a  region  where 
the  channel  narrows  to  300  ft 

•  SAM  greatly  reduces  the  sensitivity  of 
TDY  to  variations  in  conductivity  north 
of  waypoint  W.  However,  south  of  way- 
point  W,  sensitivity  is  greatly  in¬ 
creased.  Note  that  MXY  is  the  preferred 
station  triad  south  of  waypoint  W  be¬ 
cause  of  low  geometric  dilution  of 
precision 

•  SAM  has  little  impact  on  the  sensitivity 
of  TDZ  to  variations  in  conductivity. 


*TDX,  TDY  and  TDZ  are  the  respective  TDs  associated  with  the 
x,  y  and  z  secondary  stations. 
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(b)  TDY 
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Figure  2.4-3  TD  Sensitivity  to  in  (Conductivity) 

With  SAM  and  Without  SAM 
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When  these  results  are  compared  to  Fig.  2.4-2,  it  is  found 
that  Fig.  2.4-2  correctly  displays  the  effect  of  SAM  on  TD 
sensitivity  to  conductivity  at  the  waypoints.  Figure  2.4-2 
may  therefore  be  used  with  confidence,  in  the  vicinity  of 
the  river,  for  refractive  index,  the  parameter  a  and  con¬ 
ductivity. 

Although  the  analysis  herein  suggests  that  the  SAM 
can  promote  temporal  grid  instability  at  some  locations  due 
to  variations  in  propagation  parameters,  the  SAM  is  neverthe¬ 
less  critical  to  controlling  relative  errors  betweer  secondary 
and  master  station  clocks.  A  possible  remedy  of  this  situa¬ 
tion  may  be  to  establish  other  SAMs  in  the  service  area  and 
to  utilize  information  from  all  SAMs  to  derive  local  phase 
adjustments.  By  employing  optimal  estimation  theory,  it  may 
be  possible  to  reduce  grid  sensitivity  to  parameter  varia¬ 
tions  beyond  what  is  possible  with  a  single  SAM.  The  predic¬ 
tions  of  seasonal  grid  variations  presented  in  Section  3 
utilize  the  present  SAM  location  and  operating  procedure . 
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3.  PREDICTED  SEASONAL  VARIATIONS  IK 

THE  ST.  MARYS  RIVER  LORAN-C  GRID 


3.1  INTRODUCTION 

In  Section  2,  TD  sensitivity  equations  for  the  propa¬ 
gation  parameters  n,  a,  and  a  are  derived  from  theoretical 
considerations.  These  equations  permit  the  prediction  of  TD 
variations  corresponding  to  expected  propagation  parameter 
variations.  In  this  section,  historical  meteorological  data 
and  other  published  information  are  used  to  describe  the 
variation  of  each  parameter  by  a  mean  seasonal  cycle  and 
random  deviations  from  this  cycle.  Expected  mean  and  random 
TD  variations  are  then  computed,  using  the  sensitivity  equa¬ 
tions,  and  compared  with  temporal  TD  variations  observed  in 
the  St.  Marys  River  chain  data  collected  by  the  U.S.  Coast 
Guard.  The  impact  of  predicted  seasonal  TD  variations  on 
a  Loran-C  user  employing  the  calibrated  TD  grid  is  given  in 
terms  of  expected  position  errors  at  river  waypoints. 


3.2  EXPECTED  VARIATIONS  IN  PROPAGATION  PARAMETERS 

The  propagation  parameters  n,  a,  and  a  vary  in  space 
and  time  in  the  St.  Marys  River  Loran-C  service  area.  In 
order  to  satisfy  the  requirement  of  a  homogeneous  medium 
(adopted  in  Section  2  for  tractability ) ,  spatial  averages 
for  the  parameters  must  be  determined.  Secticns  3.2.1  to 
3.2.3  present  seasonal  variations  of  the  spatial  averages. 
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3.2.1  Surface  Refractive  Index 

The  refractive  index  n,  which  impacts  primary  phase 

delay  as  discussed  in  Section  2.3.1,  can  be  determined  from 

surface  meteorological  parameters.  Since  the  refractive 

index  of  air  does  not  exceed  unity  by  more  than  a  few  parts 
4 

in  10  ,  it  is  convenient  to  consider  instead  the  ref ractivity 

N  =  (n-1)  x  106  (3.2-1) 


which  is  of  the  order  of  300  units.  Ref ractivity  is  related 
to  atmospheric  pressure  p  (mbar),  absolute  temperature  T 
(°  Kelvin),  and  relative  humidity  RH  (%)  by  the  formula  (Ref.  7): 


N 


77.6 


+  3730.0 


(3.2-2) 


where  e  denotes  saturation  water  vapor  pressure  at  temperature 

b 

T.  The  first  and  second  terms  of  Eq.  3.2-2  correspond  to  the 
"dry"  and  "wet"  components  of  N,  respectively. 


Data  from  the  National  Weather  Service  Station  at 
Sault  Sainte  Marie,  Michigan  are  assumed  to  be  representative 
of  weather  conditions  in  the  St.  Marys  River  region  for  pur¬ 
poses  of  dete  i^ini  ,ig  seasonal  variations  of  ref  ractivity . 

Equation  3.2-2  has  been  employed  in  Ref.  8  to  compute  monthly 
statistics  for^N  at  Sault  Sainte  Marie  (by  hour)  from  meteoro¬ 
logical  data  collected  over  an  eight-year  period.  In  Fig.  3.2-1, 
the  mean  seasonal  cycle  and  llo  bounds  along  with  the  extreme 
bounds  are  plotted  for  12  o'clock  noon  local  standard  time. 

These  data  reveal  a  dist  inct  seasonal  cycle  with  a  relatively 
narrow  la  spread.  Furthermore,  refractivity  statistics  versus 
time  of  day  for  each  month  (August,  for  example,  is  given  in 
Fig.  3.2-2)  show  that  seasonal  variations  are  only  slightly 
dependent  on  the  time  oi  day. 
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3.2-1  Historical  Seasonal  Variations  of 
Refractivity  at  Sault  Sainte  Marie 
Weather  Station 


Figure  3.2-2 


Historical  Diurnal  Variations  of 
Refractivity  at  Sault  Sainte  Marie 
During  August 
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The  average  refractivity  for  the  18-day  September 
1977  data  collection  period  is  324.0  at  Sault  Saints  Marie, 
compared  to  the  computed  mean  of  331.0  at  the  data  collection 
sites.  Sault  Sainte  Marie  data  is  employed  consist sntly 
herein  because  continuous  weather  data  are  not  available 
at  the  sites.  Note  that  although  neither  Sault  Sainte  Marie 
nor  the  sites  necessarily  exhibit  the  average  refractivity 
of  the  coverage  area,  the  refractivity  is  within  the  la 
bounds  associated  with  the  month  of  September  shown  in  Fig. 
3.2-1. 


3.2.2  Parameter  a 


The  propagation  parameter  a,  which  affects  secondary 
phase  delay  as  discussed  in  Section  2.3.2,  depends  on  the 
vertical  lapse  rate  of  the  refractive  index,  i.e.,  dn/dh. 

The  lapse  rate  is  commonly  interpreted  as  the  change  in  n 
between  the  surface  and  a  1  km  altitude.  If  the  corresponding 
change  in  refractivity  is  AN,  then  Eq.  2.3-5  becomes 

a  =  1  +  0.006363  (3.2-3) 

For  most  meteorological  conditions  (excluding  temperature 
inversions),  AN  is  negative  and  highly  correlated  with  surface 
refractivity.  In  Ref.  9,  the  following  regression  equation 
was' fit  to  888  pairs  of  data  from  45  U.S.  weather  stations 


AN  =  -7.32  e0'005577  N 


(3.2-4) 


The  resulting  high  correlation  coefficient  (0.93)  justifies 
computing  AN  from  surface  meteorological  data  alone,  in  the 
absence  of  radiosonde  data  from  1  km  altitude. 


Although  radiosonde  data  at  Sault  Sainte  Marie  are 
employed  in  Ref.  7  to  determine  a  mean  seasonal  cycle  for 
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AN,  the  data  records  necessary  to  determine  lo  and  extreme 
bounds  are  not  readily  obtainable.  Therefore,  Eq.  3.2-4  is 
employed  to  transform  the  statistics  of  surface  refractivity 
to  statistics  for  AN,  resulting  in  the  seasonal  variations 

t 

plotted  in  Fig.  3.2-3.  The  mean  seasonal  cycle  computed  in 
this  manner  agrees  closely  with  the  radiosonde  data  of  Ref.  7 
The  average  value  of  AN  at  Sault  Sainte  Marie  for  the  Sep¬ 
tember  1977  data  collection  period  is  -45. 


R-37«3« 


MONTH 


Figure  3.2-3  Seasonal  Variations  of  AN 

at  Sault  Sainte  Marie 

3.2.3  Ground  Conductivity 

Unlike  the  atmospheric  parameters  n  and  a,  ground 
conductivity  ac  (mho/m)  cannot  be  easily  related  to  physical 
parameters  and  must  be  inferred  from  signal  propagation  data. 
Although  no  known  published  conductivity  data  is  available 
near  Loran-C  frequencies  (90-110  KHz)  for  the  St.  Marys  River 
region,  data  at  other  frequencies  are  reported  in  Refs  10-13 


THE  ANALYTIC  SCIENCES  CORPORATION 


Since  conductivity  varies  little  with  frequency  below  10  MHz 
(Ref.  10),  published  data  are  employed  to  establish  an  expected 
conductivity  range  of  0.001  mho/m  to  0.01  mho/m  for  the  St. 
Marys  River  Loran-C  application.  The  "average"  ground  conduc¬ 
tivity  of  interest  here  is  that  associated  with  station-to- 
user  signal  propagation  paths.  Since  determination  of  the 
average  conductivity  for  the  region  from  the  limited  infor¬ 
mation  cited  above  is  impractical  and  unwarranted,  the  conduc¬ 
tivity  has  been  estimated  from  the  September  1977  Loran-C  data. 

The  parameters  n  and  a  are  assumed  constant  at  their 

known  average  values  for  the  period,  and  a  is  estimated  using 

the  homogeneous  Loran-C  signal  propagation  model  established 

in  Section  2.3.  A  value  of  0.0047  mho/m  for  the  parameter 

o  has  been  identified  from  the  September  1977  U.S.  Coast 
c 

Guard  TD  data  by  maximum  likelihood  estimation  techniques 
(Ref.  14)  which  are  conveniently  embodied  in  TASC  software 
PARA IDE.  The  resulting  homogeneous  model  is  not  a  replace¬ 
ment  for  the  RB  model,  as  it  yields  much  larger  residual  TDs. 
However,  the  calibrated  conductivity  (0.0047  mho/m)  serves  as 
a  baseline ,  to  be  used  with  other  published  information  to 
provide  expected  seasonal  conductivity  variations.  It  must 
be  remembered  that  0.0047  mho/m  has  been  derived  from  data 
and,  therefore,  is  "tuned"  to  the  specific  propagation  con¬ 
ditions  inherent  in  the  data.  In  addition,  the  estimate 
may ^ be  influenced  by  measurement  noise. 

Since  the  spatial-average  conductivity  for  the  region 
is  not  expected  to  vary  appreciably  except  in  winter,  the 
mean  of  the  spatial-average  conductivities  over  man/  years 
is  assumed  equal  to  0.0047  mho/m  for  spring,  summer  and  fall. 
The  random  variability  in  the  conductivity  from  yea:1  to  year 
for  these  three  seasons  is  estimated  by  assuming  that  the 
published  range  of  conductivities  over  space  (0.001  to  0.01 
mho/m)  is  the  same  range  encountered  in  the  spatial- average 
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• 

conductivity  from  year  to  year.  If  this  range  defines  a 
3a  spread  for  spring,  summer  and  fall  conductivities,  then 
the  la  factor*  for  conductivity  is  1.5  and  the  la  spread 
is  0.0031  to  0.0071  mho/m. 

A  significant  decrease  in  conductivity  is  expected 

during  winter  in  the  St.  Marys  River  region  due  to  snow  and 

-5  -4 

ice  («  3.0  x  10  mho/m)  and  frozen  ground  (w  10  mho/m), 

which  have  lower  conductivities  (Ref.  15)  than  average  ground 

(a  0.001  mho/m).  A  two-layer  earth  model  (Ref.  16)  has  been 

employed  to  estimate  the  effective  conductivity  for  1-2  m 

of  snow  or  frozen  ground  on  top  of  average  ground,  and  1  m 

of  ice  on  fresh  water.  The  minimum  calculated  conductivity 

was  0.00057  mho/m,  for  snow  on  average  ground.  Thus,  the  range 

of  conductivities  over  the  region  for  all  winters,  including 

mild  winters,  is  expected  to  be  0.00057  to  0.01  mho/m.  The 

geometrical  mean  of  the  extremes,  0.0024  mho/m,  is  employed 

as  the  mean  spatial-average  conductivity  for  winter.  Since 

the  extent  of  ice,  snow  and  frozen  ground  can  vary  greatly 

from  winter  to  winter,  the  expected  range  is  interpreted  only 

as  a  2o  spread.  The  resulting  la  conductivity  factor  is  2 

and  the  la  spread  is  0.0012  to  0.0048  mho/m.  The  seasonal 

variations  in  the  TD  grid  presented  in  the  next  section  should 

be  recomputed  if  better  estimates  of  conductivity  variations 

can  be  determined. 

3.3  EXPECTED  VARIATIONS  IN  TDs  * 

In  Table  3.3-1,  statistics  of  seasonal  parameter 
variations  estimated  in  Section  3.2  are  summarized  for  four 
times  of  the  year,  corresponding  approximately  to  the  middle 
of  the  four  seasons.  The  mean  cycle  is  an  a  priori  estimate 
of  seasonal  parameter  changes  while  the  standard  deviation  (la) 

*A  la  factor  of  1.5  implies  that  the  standard  deviation  of 
£n  ac  is  In  1.5.  The  la  spread  about  the  mean  oc  is  cW 
1.5  to  1.5  oc. 
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TABLE  3.3-1 

EXPECTED  SEASONAL  VARIATION  OF  PARAMETERS 


Season 

D 

a 

oc  (mho/ i) 

Mean 

la( xlO6 ) 

Mean 

lo 

Mean 

lo  Factor 

Spring 

1.000305 

11 

0.76 

0.02 

0.0047 

Summer 

1.000328 

16 

wm 

0.03 

0.0047 

mm 

Fall 

1.000312 

13 

0.73 

0.03 

0.0047 

H 

Winter 

1 . 000304 

4 

0.77 

0.01 

0.0024 

HI 

September 

1977 

1 . 000324 

- 

0.70 

- 

0.0047 

- 

indicates  the  random  fluctuation  about  the  mean  cycle 

over  many  years.  Also  indicated  in  Table  3.3-1  are  the 

parameter  values  for  the  period  during  September  1977  when  the 

TD  data  used  to  calibrate  the  RB  model  were  collected.  The 

primary  objective  is  to  present  statistics  for  TD  variations 

at  the  river  waypoints  corresponding  to  expected  variations 

in  the  propagation  parameters  n,  a,  and  a  .  The  sensitivity 

c 

equations,  which  relate  TD  variations  to  parameter  varia¬ 
tions,  were  derived  in  Section  2  and  are  summarized  in 
Appendix  A. 

3.3.1  Mean  Seasonal  TD  Shifts 

|  The  mean  TD  shift  for  each  season  is  the  expected 

change  in  TD  experienced  if  all  three  propagation  parameters 
change  from  their  values  during  the  September  1977  data  col¬ 
lection  period  to  the  assumed  mean  values  for  the  season. 

Mean  winter  TD  shifts  are  illustrated  for  downbound  river 
waypoints  in  Fig.  3.3-1,  and  mean  TD  shifts  for  all  seasons 
and  all  waypoints  are  tabulated  in  Appendix  C  (Tables  C-l  to 
C-4 ) . 
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DISTANCE  ALONG  RIVER  (mi) 

Figure  3.3-1  Predicted  Mean  TD  Shifts  in  Winter 

at  Downbound  Waypoints 


The  results  indicate  that  mean  TD  shifts  do  not 
exceed  5  nsec  for  spring,  summer  and  fall,  but  can  exceed 
200  nsec  (with  an  average  of  61  nsec)  for  winter.  This 
behavior  is  attributed  to  the  following: 


•  Phase  delay  sensitivity  to  n  (0.001  nsec/ 
ysec  per  change  of  10-6  in  n)  ard  to  a 
(0.005  nsec/ysec  per  change  of  0.01  in 

a)  is  insignificant  compared  with  sen¬ 
sitivity  to  cc. 

•  The  relative  change  in  phase  delay  for 
a  change  in  conductivity  from  0.0047 
to  0.0024  mho/m  is  approximately  0.7 
nsec/ysec. 

•  Mean  conductivities  for  spring,  summer 
and  fall  have  been  assumed  to  ecual 

the  estimated  conductivity  for  September 
1977  data  collection.  Therfore,  only 
changes  in  n  and  a  contribute  to  mean 
TD  shifts  for  these  seasons. 
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•  Resulting  mean  winter  TD  shifts  are  a 
direct  consequence  of  the  assumed  de¬ 
crease  in  mean  conductivity  between 
September  1977  and  winter. 

Although  the  specific  mean  TD  shifts  predicted  by  this  study 
can  only  be  as  reliable  as  the  assumed  mean  parameter  varia¬ 
tions,  two  conclusions  can  be  drawn  with  some  confidence. 

First ,  changes  in  refractive  index  and  the  parameter  a  do  not 
contribute  more  than  20  nsec  to  TD  variations.  Second,  there 
is  a  potential  for  changes  in  ground  conductivity  to  produce 
significant  TD  variations. 

An  additional  observation  can  be  made  upon  inspection 
of  Fig.  3.3-1.  There  is  no  mean  shift  in  TDs  at  tie  SAM 
(located  near  waypoint  L) ,  but  the  magnitude  of  the  shift 
increases  as  the  differential  range  from  SAM  (magnitude  of 
the  difference  between  waypoint  TD  and  SAM  TD)  increases. 
Although  the  shift  does  not  increase  monotonically  with  dis¬ 
tance  along  the  river  from  SAM,  the  greatest  shift  tends  to 
be  near  the  extreme  waypoints. 

3.3.2  Random  Seasonal  TD  Variations 

In  any  particular  year,  the  actual  seasonal  cycle  of 
a  propagation  parameter  may  differ  significantly  from  the 
mean  cycle,  which  is  the  average  over  many  years.  The  year- 
to-year  parameter  variations  are  referred  to  here  as  random 
variations  because  they  cannot  be  identified  a  pric  r i .  There 
may  also  be  variations  within  a  season  in  a  particular  year, 
but  the  available  data  do  not  permit  the  statistics  of  these 
variations  to  be  determined.  Expected  random  (lo)  TD  varia¬ 
tions  corresponding  to  the  assumed  random  (lo)  part  meter 
variations  of  Table  3.3-1,  are  illustrated  for  sumr er  and 
winter  in  Figs.  3.3-2  and  3.3-3.  Random  TD  variations  for 
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Figure  3.3-1:  Predicted  Random  TD  Variations  in 

Summer  at  Downbound  Waypoints 


Figure  3.3-3  Predicted  Random  TD  Variations  in 

Winter  at  Downbound  Waypoints 

all  seasons  and  all  waypoints  are  tabulated  in  Appendix  C 
(Tables  C-l  to  C-^ ) . 

Because  \ariations  in  n  and  a  contribute  little  to 
TD  variations  (see  Section  3.3.1)  and  because  random  conduc¬ 
tivity  variations  for  spring,  summer  and  fall  are  assumed 
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equal,  random  TD  variations  are  essentially  the  sama  for 
these  seasons  (approximately  125  nsec  maximum  and  50  nsec 
rms).  Random  TD  variations  for  winter  increase  to  approxi¬ 
mately  200  nsec  maximum  and  75  nsec  rms,  due  to  the  assumed 
increase  in  the  lo  conductivity  factor  (defined  in  Section 
3.2.3)  from  1.5  in  other  seasons  to  2.0  in  winter.  On  compar¬ 
ing  Figs.  3.3-1  and  3.3-3,  it  is  observed  that,  in  winter, 
the  mean  TD  shifts  approximately  equal  the  random  TD  varia¬ 
tions.  If  a  priori  compensation  for  mean  TD  shifts  were 
applied  to  the  calibrated  grid,  then  random  temporal  grid 
instability  would  remain.  If  the  random  component  of  grid 
instability  varys  slowly  over  the  year,  then  compensation 
may  be  determined  and  applied  periodically.  If  the  variations 
are  rapid,  however,  then  some  form  of  real-time  compensation 
would  be  required.  The  results  presented  in  this  section  are 
based  on  reasonable  estimates  of  conductivity  variations  and 
are  subject  to  revision  if  better  conductivity  data  becomes 
available . 

3.3.3  Comparison  With  Observed  Variations 

Subsequent  to  collecting  data  during  September  1977 
for  grid  calibration,  the  U.S.  Coast  Guard  revisited  10  of 
the  data  collection  sites  as  many  as  four  times  between 
November  1977  and  March  1978.  Figure  3.3-4  shows  site  TDs 
(after  editing  outliers  and  averaging  an  ensemble  of  0.5 
hr  data  collection  intervals)  relative  to  TDs  predicted  by 
the  RB  model .  for  the  10  revisited  sites  (see  Fig.  3.3-5) 
and  all  five  data  collection  periods.  This  chart  reveals 
temporal  variations  in  site  TDs,  as  well  as  differences 
(residuals)  between  observed  TDs  and  the  calibratec  grid. 

The  residuals  for  September  1977  (data  collection 
period  A)  may  include  unmodelled  spatial  details  (local 
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Figure  3.3-4 


TD  Data  Relative  to  Calibrated  Grid 
for  Five  Data  Collection  Periods 
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Figure  3.3-5  Data  Collection  Sites 

warpages),  temporal  grid  instability  over  the  18  davs  of  data 
collection,  and  measurement  errors.  Temporal  variations  over 
18  days  are  expected  to  be  small  based  on  observed  variations 
in  refractive  index  and  the  parameter  a  and  on  reasonable 
variations  in  conductivity.  (However,  some  temporal,  variation 
may  be  caused  by  weather  effects  such  as  rain  or  th?  passage  of 
weather  fronts.)  Significant  local  warpages  appareitly  exist, 
as  indicated  by  nearly  constant  residuals  for  all  dvta  col¬ 
lection  periods  at  some  sites,  such  as  113  (TDX)  and  108  (TDZ). 
The  relative  importance  of  temporal  variations,  loc  il  warpages 
and  measurement  errors  cannot  be  determined  from  the  present 
data  .set  .  It  is  conceivable  that  calibration  residuals  would 
be  smaller  than  those  shown  in  Fig.  3.3-4  if  the  da. a  did  not 
include  temporal  variations  and  measurement  errors. 
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The  TD  variations  between  data  collection  periods 
revealed  by  Fig.  3.3-4  are  of  particular  interest  to  this 
study.  Site  TDs  corresponding  to  the  data  collection  per¬ 
iods  cannot  be  predicted  without  first  relating  conductivity 
to  physical  ground  parameters  for  the  periods.  However,  a 
comparison  of  observed  variations  to  those  predicted  in 
Sections  3.3.1  and  3.3.2  reveals  the  following: 

•  SAM  TDs  (site  124)  vary  less  thar;  10  nsec 

•  TD  variations  at  sites  102,  123,  and 
113  are  under  100  nsec  peak-to-pe ak 
and  may  be  explainable  by  expected 
parameter  variations 

•  TD  variations  at  the  other  six  sites 
are  unusually  large  (850  nsec  ma> i- 
mum  and  300  nsec  average)  and  carnot 
be  easily  explained  by  assumed  param¬ 
eter  variations 

•  The  average  peak-to-peak  TD  varirtion 
over  all  revisited  sites  is  240  i sec . * 

The  unexplained  temporal  variations  may  be  caused  by  some¬ 
thing  other  than  propagation  effects  such  as  d:  fferences  in 
the  receiving  antenna  ground  plane  during  the  separate  data 
collection  periods,  or  unmodelled  weather  effects.  It  could 
also  be  that  some  sites  such  as  those  identified  as  anomalous 
in  Ref.  1  (namely  104,  108,  and  117)  may  exhibit  localized 
temporal  variations.  The  large  variations  at  sites  118  (X 
station)  and  119  (Y  station)  were  very  likely  caused  by  a  change 
in  the  grounding  environment  associated  with  moving  the  re¬ 
ceiving  antenna  from  the  top  of  the  transmitting  tower  (in 
September)  to  the  earth  (in  November).  Although  observed  TD 
variations  apparently  cannot  be  explained  by  enpected  variations 


^Editing  of  selected  data  at  sites  104,  117,  1:8,  119  and  121, 
which  may  be  anomalous,  reduces  the  average  tc  100  nsec. 
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in  the  propagation  parameters  considered  herein,  this  should 
be  substantiated  by  a  long-term  signal  monitoring  effort. 


3.4  POSITION  ERRORS  AT  WAYPOINTS 

In  Section  3.3,  seasonal  statistics  were  determined 
for  the  component  of  TD  error  caused  by  changes  in  the  propa¬ 
gation  parameters.  These  results  are  employed  in  this  section 
to  predict  user  position  errors  at  river  waypoints  using  the 
transformation  equations  given  in  Appendix  B. 

If  the  calibrated  grid  is  not  compensated  for  mean 
winter  TD  shifts  (see  Section  3.3.1),  then  a  Loran-C  user 
employing  the  grid  in  winter  is  expected  to  experience  a 
deterministic  (mean)  position  offset.  The  magnitude  and  bear¬ 
ing  of  the  mean  position  errors  at  the  waypoints  are  tabu¬ 
lated  in  Appendix  C  (Table  C-5)  for  the  triad  of  transmitter 
stations  which  yields  minimum  error  with  the  calibrated  grid 
model.  Although  the  magnitude  of  the  mean  winter  offset  has 
an  average  value  of  approximately  60  ft  (over  all  waypoints), 
a  correction  table  could  be  employed  by  the  user  to  compensate 
for  this  offset  in  the  winter  TD  grid. 

In  addition  to  mean  TD  shifts,  random  seasonal  grid 
instability  (see  Section  3.3.2)  also  contributes  to  position 
err  ir.  This  contribution  is  plotted  in  Fig.  3.4-1  for  summer* 
and  for  winter  at  the  downbound  waypoints  for  the  same  station 
triad  as  used  to  evaluate  the  calibrated  grid  in  Ref.  1  (see 
Table  C-5).  Random  position  errors  are  tabulated  in  Appendix 
C  (Table  C-6)  for  all  station  triads.  This  random  component 


^Random  propagation  errors  for  spring  and  fall  are  .essentially 
the  same  as  for  summer. 
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H  38013 


Figure  3.4-1  Component  of  2d  rms  Position  Error 

at  Downbound  Waypoints  Due  to 
Random  Seasonal  Grid  Instability 


of  grid  instability  is  interpreted  as  the  variability  in 
the  grid  from  one  year  to  another  during  the  specified  season. 
The  variability  over  a  given  season  is  difficult  to  assess 
since  historical  meteorological  data  are  only  available  as 
monthly  averages  and  conductivity  data  are  no'  available. 

The  increase  in  position  error  at  the  extremes  of  the  river 
is  due  to  the  combined  effect  of  SAM  and  GDOP*.  Note  that 
even  if  the  RB  model  provides  an  error-free  prediction  of 
the  grid  which  existed  during  September  1977  data  collection, 
there  would  still  be  2d  rms  position  errors  (Eq.  B-6)  of  75  ft 
(rms  over  all  waypoints)  in  spring,  summer  and  fall  and  130  ft 

^Geometric  Dilution  of  Precision 
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winter  from  one  year  to  another.  However,  it  should  also  be 
noted  that  this  predicted  variation  in  the  grid  is  due  primar¬ 
ily  to  predicted  variations  in  conductivity  from  year  to  year. 
The  practicality  of  compensating  (by  the  user)  for  these 
variations  in  the  grid  depends  on  the  frequency  of  the  fluctua¬ 
tions.  For  example,  if  the  variation  is  slow,  it  may  be 
possible  to  provide  a  simple  first-order  grid  compensation 
algorithm  in  the  user  equipment  which  is  valid  for  a  month 
or  even  the  entire  season.  The  result  would  be  to  reduce  the 
predicted  errors  due  to  the  modelled  random  effects. 

In  actuality,  the  RB  model  does  not  predict  the 
September  1977  grid  exactly  and  other  contributors  to  posi¬ 
tion  error  must  be  considered.  Equation  2.2-2  expresses  the 
total  TD  error  invoked  when  using  the  RB  model  to  predict 
the  true  TD  (the  grid  measured  in  the  absence  of  noise).  The 
total  error  includes  three  components: 

•  Residual  grid  prediction  error  (model 
calibration  error)  —  error  due  to  an 
inability  to  determine  exactly  the 
coefficients  of  the  RB  model  with  a 
finite  amount  of  noisy  data 

•  Local  warpage  —  unmodelled  spatial 
details  which  do  not  change  with  time 

•  Temporal  grid  instability  —  varia¬ 
tions  due  to  changing  propagation 
conditions  (seasonal,  diurnal,  and 
weather  front  passage  effects). 

If  position  is  established  at  a  waypoint  by  employing  noise¬ 
less  TD  measurements  and  the  calibrated  RB  grid,  then  the 
total  position  error  (relative  to  the  actual  waypoint)  in¬ 
cludes  a  component  due  to  each  of  the  above  components  of 
TD  error.  In  Fig.  3.4-2  the  residual  grid  prediction  (position) 
error  is  reproduced  from  Ref.  1.  This  curve  is  a  theoretical 
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lower  bound  on  the  expected  2d  rms  position  error  and  can  only 
be  approached  if  both  local  warpage  and  temporal  grid  insta¬ 
bility  are  negligible.  It  should  be  kept  in  mind  that  this 
bound  is  a  characteristic  of  the  particular  spatial-area 
model  structure  chosen  (range  and  bearing  polynomial)  and 
of  the  quality  and  quantity  of  the  particular  data  used  for 
calibration. 


R  3*012 


Figure  3  4-2  Residual  Grid  Prediction  Errors  and 

Total  September  1977  and  Winter 
Errors  at  Downbound  Waypoints 


An  upper  bound  on  the  expected  2d  rms  position  error 
at  the  waypoints  for  the  September  1977  data  collection  period 
is  shown  in  Ref.  1  and  reproduced  in  Fig.  3.4-2.  This  bound 
was  arrived  at  on  an  ad  hoc  basis  in  Ref.  1  in  an  attempt  to 
characterize  the  observed  site-to-site  variability  of  the 
September  1977  data  relative  to  the  RB  model.  Local  warpage 
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and  temporal  grid  instability,  as  well  as  measurement  uncer¬ 
tainties  (e.g.,  potential  site-to-site  ground  plane  differences), 
all  contribute  to  the  observed  site-to-site  variability  in  the 
data.  Additional  data  collected  since  the  September  1977  data 
collection  period  are  not  sufficient  to  allow  this  bound  to 
be  changed  at  this  time.  A  Loran-C  user,  employing  the  cali¬ 
brated  grid  on  the  St.  Marys  River  during  the  fall  of  another 
year,  will  still  encounter  local  warpage  due  to  any  unmodelled 
spatial  details  of  signal  propagation.  However,  the  user  would 
probably  not  experience  the  same  measurement  problems  or  short¬ 
term  temporal  effects  which  may  have  occurred  during  September 
1977  data  collection.  On  the  other  hand,  the  contribution  of 
fall  grid  instability  (from  year  to  year,  as  discussed  in  this 
report)  to  user  position  error  is  expected  to  be  greater  than 
the  contribution  of  grid  instability  over  the  18  days  of 
September  1977  data  collection.  Because  the  relative  con¬ 
tribution  of  local  warpage,  temporal  grid  instability  and 
measurement  error  to  the  predicted  September  1977  position 
errors  is  not  known,  the  total  user  position  error  during 
another  fall  cannot  be  determined  exactly. 

If  the  predicted  position  errors  for  September  1977 
data  collection  are  assumed  to  be  representative  of  errors 
for  any  fall,  however,  then  these  errors  can  serve  as  a  base¬ 
line  for  predicting  winter  position  errors.  The  total  winter 
position  errors  shown  in  Fig.  3.4-2  (and  tabulated  in  Table 
C-7)  were  determined  by  employing  this  fall  baseline  along 
with  the  incremental  change  in  predicted  position  errors  be¬ 
tween  fall  and  winter  due  to  random  grid  instability  (shown 
in  Fig.  3.4-1).  The  total  winter  position  errors  should  be 
interpreted  as  an  expected  upper  bound  on  the  predicted  posi¬ 
tion  error  for  a  Loran-C  user  navigating  the  St.  Marys  River 
in  any  season.  It  is  evident  from  Fig.  3.4-2  that  the  addi¬ 
tional  contribution  of  winter  grid  instability  to  position 
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error  is  small  when  combined  in  a  root-sum-square  sense  with 
position  errors  for  the  September  1977  data  collection  period 
predicted  in  Ref.  1.  The  2d  rms  position  error  (rms  over  all 
waypoints)  increases  from  290  ft  for  September  1977  data  col¬ 
lection  to  310  ft  for  winter. 

When  interpreting  the  total  position  errors  predicted 
herein,  it  must  be  remembered  that  the  observed  site-to-site 
variability  in  the  1977  data,  relative  to  the  RB  model,  was 
the  basis  for  the  mathematical  model  used  to  predict  the  asso¬ 
ciated  total  position  error  at  the  waypoints.  If  a  major  com¬ 
ponent  of  this  observed  site-to-site  variability  is  due  to 
local  warpage  effects  which  are  constant  with  time,  it  may  be 
possible  to  reduce  these  errors  (e.g.,  by  using  force-fit 
techniques)  and  thereby  reduce  the  total  error  at  the  waypoint 
Because  the  resulting  total  error  would  be  smaller,  the  rela¬ 
tive  change  between  fall  and  winter  would  become  more  signifi¬ 
cant.  A  comprehensive  data  monitoring  effort  should  permit 
the  major  contributors  to  total  grid  error  to  be  isolated  and 
thereby  allow  better  predictions  of  position  error  to  be 
determined. 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 


4.1  CONCLUSIONS 

This  study  has  predicted  the  impact  of  expected  sea¬ 
sonal  variations  in  propagation  parameters  on  the  stability  of 
the  St.  Marys  River  Loran-C  time  difference  grid.  Sensitivity 
of  the  grid  to  the  important  propagation  parameters  —  surface 
refractive  index  n,  the  parameter  a  associated  with  refractive 
index  gradient,  and  ground  conductivity  oc  —  follows  from 
classical  propagation  theory.  TD  variations,  corresponding  to 
estimates  of  seasonal  parameter  variations,  were  computed 
and  compared  with  TD  variations  embodied  in  U.S.  Coast  Guard 
data.  Position  errors  at  river  waypoints  were  then  determined 
from  the  TD  variations  in  order  to  assess  the  impact  of  tem¬ 
poral  grid  instability  on  a  Loran-C  user. 

During  the  study,  it  was  found  that  seasonal  variations 
in  the  parameters  n  and  a  can  be  characterized  statistically 
from  historical  meteorological  data  recorded  at  the  National 
Weather  Service  Station  at  Sault  Sainte  Marie,  Michigan. 

However,  characterization  of  the  ground  conductivity  is  compli¬ 
cated  by  the  heterogeneity  of  the  region  and  by  the  inadequacy 
of  published  data.  The  following  seasonal  analysis  results 
are  based  on  best  available  estimates  of  seasonal  parameter 
variations : 

•  The  effect  of  variations  in  n  and  a  on 
grid  instability  is  less  than  20  nsec 
at  all  waypoints 

•  The  mean  difference  between  TDs  in  spring, 
summer,  and  fall,  and  TDs  in  winter  is 

60  nsec  averaged  over  all  waypoints 
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•  Random  (la)  TD  variations  (about  the 
mean)  due  to  assumed  conductivity 
variations  are  approximately  45  nsec 
rms  for  spring,  summer  and  fall  and 
75  nsec  rms  for  winter. 


The  predictions  of  grid  instability  computed  in  this  study 
cannot  explain  the  TD  variations  observed  at  most  data  col¬ 
lection  sites  (240  nsec  peak-to-peak,  averaged  over  all  re¬ 
visited  sites) .  This  conclusion  is  based  upon  theoretical 
models  and  expected  propagation  parameter  variations,  and 
should  be  validated  by  a  more  comprehensive  data  collection 
effort . 


Predicted  position  errors  at  the  waypoints  have  been 
computed  corresponding  to  both  mean  seasonal  TD  shifts  and 
random  (la)  TD  variations.  The  deterministic  (mean)  position 
offsets  corresponding  to  mean  winter  TD  shifts  average  60  ft 
in  magnitude  but  can  be  removed  by  grid  compensation.  Random 
position  errors  (Fig.  3.4-1)  cannot  be  removed  by  a  priori 
grid  compensation,  and  add  to  position  errors  due  to  grid 
prediction  (calibration)  and  local  warpage  (unmodelled  spatial 
details).  If  the  RB  model  were  capable  of  perfectly  repre¬ 
senting  the  grid  which  existed  during  September  1977,  then 
estimated  random  seasonal  grid  instability  would  still  contribute 
an  average  (over  all  waypoints)  of  75  ft  2d  rms  position  error 
in  spring,  summer  and  fall  and  130  ft  in  winter.  An  average 
total  2d  rms  position  error  of  290  ft  associated  with  the 
upper  bound  on  position  error  during  September  1977  data  col¬ 
lection  (Ref.  1  and  Fig.  3.4-2),  and  assumed  to  be  representa¬ 
tive  of  conditions  occurring  during  any  fall,  increases  to 
310  ft  for  winter.  It  should  be  kept  in  mind,  however,  that 
the  predicted  waypoint  position  errors  in  Ref.  1  may  include 
the  effect  of  site-to-site  measurement  errors  (e.g.,  antenna 
ground  plane  differences)  which  cannot  be  identified  from  the 
available  data  and  would  not  be  experienced  by  a  user  at  the 
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waypoints.  Also,  there  is  some  evidence  from  the  revisited 
site  data  that  local  warpage  may  be  a  significant  contributor 
to  the  observed  site-to-site  grid  prediction  errors.  Cali¬ 
bration  of  these  local  warpages  (assuming  they  do  not  change 
with  time)  would  also  reduce  the  predicted  waypoint  position 
errors.  Finally,  ground  conductivity  variations  constitute 
the  major  contribution  to  predicted  seasonal  grid  instability. 
If  these  variations  in  conductivity  are  relatively  slow,  it 
may  be  possible  to  provide  compensation  in  the  user  equipment 
which  is  valid  for  a  month  or  an  entire  season . 


4 . 2  RECOMMENDATIONS 

In  order  to  quantify  the  apparent  grid  instability 
observed  in  the  St.  Marys  River  Loran-C  data,  it  is  recom¬ 
mended  that  a  comprehensive  data  collection  el  fort  be  ini¬ 
tiated.  Preliminary  considerations  suggest  the  following 
approach: 

•  Isolate  receiver/ environment  instabilities 

•  Collect  long-term  data  at  strategically- 
located  sites 

•  Continue  collecting  short-term  seasonal 
data  at  other  translocation  sites 

•  Analyze  long-term  and  short-tern  data  to 
identify  grid  instability. 

Each  of  the  above  is  expanded  upon  below. 

It  is  first  necessary  to  determine  wl  ether  or  not  the 
data  collection  procedure  itself  has  contributed  to  the  appar¬ 
ent  instability.  The  adequacy  of  the  receivirg  antenna  ground 
plane  should  be  determined,  especially  as  influenced  by  site 
environment,  e.g.,  1-2  m  deep  snow  cover.  Furthermore,  the 
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potential  impact  of  precipitation  and  fog  (at  the  antenna)  on 
the  phase  of  the  received  signal  should  be  tested. 

It  is  recommended  that  long-term  time  of  arrival  (TOA) 
data  be  collected  continuously  (%  15-minute  averages)  for  sev¬ 
eral  months  at  a  few  strategically-located  sites,  to  be  deter¬ 
mined  by  covariance  simulation.  TOA  data  is  especially  useful 
in  isolating  those  propagation  effects  which  are  not  experienced 
by  the  entire  chain.  In  addition  to  continuous  long-term  data, 
it  is  suggested  that  four-hour  TD  data  be  collected,  at  least 
seasonally,  at  translocation  sites  where  long-term  data  are 
not  collected.  These  data,  in  combination  with  long-term 
data,  would  permit  the  determination  of  the  spatia]  extent 
of  any  observed  instability. 

In  order  to  interpret  the  data,  the  following  additional 
information  would  be  valuable: 

•  SAM-applied  local  phase  adjustments 

•  Signal  to  noise  ratio,  ECD,  and  re¬ 
ceiver  settings 

•  Meteorological  data. 

Analysis  of  data  from  m  intensive  collection  program,  such  as 
the  above,  should  reveal  the  extent  of  temporal  grid  instability 
and  its  relative  impact  on  the  operational  status  of  the  St. 

Marys  River  Loran-C  chain.  If  the  apparent  grid  instability 
is  substantiated  by  this  data  analysis,  the  data  should  be  used 
to  establish  the  feasibility  of  compensation  algorithms  for 
temporal  grid  instability. 
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APPENDIX  A 

COMPUTATION  OF  TD  VARIATIONS 


The  sensitivity  equations  developed  in  Sections  2.3 
and  2.4  can  be  employed  to  compute  mean  TD  shifts  (Section 
3.3.1)  and  random  (la)  TD  variations  (Section  3.3.2)  at  any 
location  in  the  St.  Marys  River  Loran-C  signal  coverage  area. 
These  equations  are  summarized  below. 

Mean  TD  shifts  follow  from  Eqs.  2.3-2,  2.3-6,  2.3-7, 
2.4-2,  and  2.4-3.  The  mean  shift  (expected  value,  £)  in  TD.^ 
(usee)  from  September ' 1977  (calibration)  to  another  time 
(season)  is 

£UTD.)  *  £(TDt)  season  -  SCTD,)  calibratlon  (A- 

where  mean  TDs  are  computed  from 


tfCTOj)  -  [(W  ♦  (SF1-SFal)]  -  |(Ti-T->  ♦  (SFi-SF;))  *  TD° 


(A-2) 


Note  that  TD° ,  the  reference  TD  for  SAM,  does  not  influence 
mean  TD  shifts.  In  Eq .  A-2,  primary  and  secondary  phase  delays 
are  denoted  by  T  and  SF .  respectively;  unprimed  phase  delays 
are  between  the  transmitting  stations  and  the  user,  while  primed 
(')  phase  delays  are  between  the  stations  and  SAM;  and  the  sec¬ 
ondary  and  master  stations  are  distinguished  by  subscripts  i 
and  m.  Phase  delays  are  computed  from  mean  propagation  parameters 
n,  a,  and  c  (for  calibration  or  season)  by 


T 


R 


( usee ) 


(A-3) 
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2  1 

SF  =  I  I  A 

ks— 2  £=— i 


kH  (ln  ac)  T  +  0-0005  (a-0.75)  T 


( A-4 ) 


where 


c 

A 

R 


k£ 


983.5690892  ft/ysec 
coefficients  from  Table  2.3-1 

path  distance  (ft)  from  a  station  to 
the  user  or  SAM. 


Note  that  Eq.  A-4  is  only  valid  for  mean  conductivities  0.001 
to  0.01  mho/m. 


Random  (lo)  TD  variations  follow  from  Eqs.  2.4-2,  2.4-4 


2.4-5,  and  2.4-6.  The  la  TD  variation,  ATD.^ ,  corresponding  to 


la  variations  of  An,  Aa ,  and  A^  (where  5  *  £n  c„ )  is 
ATD.  =  [(h*  An)2  +  (h£  Acx  )2  +  (h*  A^)2  ] 


c ' 
1/2 


(A-5) 


where  the  required  sensitivity  coefficients  are  defined  below: 


hn  "  <VV  -  (Ti-Tm> 


(  A-6 ) 


h^  =  0.0005  [(T.-Tm)  -  (T^-T^)J 


( A- 7 ) 


hi  *  -0.039  f ( T . a-T  a)  -  (T:a-T,a)l 
l  l  m  i  m  '  \ 


0.47 


( A-8 ) 
(  A-9 ) 


Note  that  Eq .  A-5  permits  the  contribution  of  each  parameter 
variation  to  be  isolated. 
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APPENDIX  B 

COMPUTATION  OF  POSITION  ERRORS 


Equations  for  transforming  TD  errors  to  position 
errors,  similar  to  equations  developed  in  Ref.  1,  are  pre¬ 
sented  below.  If  a  hyperbolic  position  fix  is  obtained  using 
secondary  stations  i  and  j ,  then  the  relationship  between  TD 
errors  (ATD.^  and  ATD ^ )  and  north  (AN)  and  east  (AE)  position 
errors  is 


AN 

=  H 

ATD . 

l 

AE 

u 

ATD  . 

L  J  J 

where  Hu  is  a  2x2  matrix  which  depends  on  station  and  user 
locations  and  is  related  to  the  crossing  angle  between  the 
hyperbolic  lines  of  position.  The  mean  (expected  value,  g 
position  error  corresponding  to  mean  TD  shifts  (see  Appendix 
A,  Eq.  A-l)  is  thus  given  by 


e  (an)' 

=  H„ 

<?  (ATDi) 

£  (AE) 

u 

€  (ATD  ) 

►  J  ■ 

and  may  be  interpreted  as  a  vector  with  magnitude  and  bearing. 


In  addition  to  mean  TD  shifts,  there  are  random 
position  errors  due  to  random  variations  in  propagation  para¬ 
meters  as  well  as  due  to  residual  grid  prediction  error,  local 
warpage ,  and  user  measurement  noise.  Random  position  errors 
are  represented  by  the  covariance  matrix  of  the  position  error 
vector  (AN,  AE).  The  total  position  error  covariance  for 
zero  measurement  noise  is  computed  (see  Ref.  17)  by 


B-l 


A: 
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P  =  H  [R  +  R  +  R  1  HT 
u  L  rg  iw  tgj  u 

where  the  superscript  T  denotes  transpose  and 


( B-3 ) 


R  =  covariance  of  residual  grid  prediction 
g  error 

R^w  =  covariance  of  local  warpage  erro:’ 

Rt  =  covariance  of  temporal  grid  instability  . 

Rrg  is  related  to  the  error  covariance  of  the  RB  model  coeffi¬ 
cients,  and  depends  only  on  the  quantity  and  quality  of  cali¬ 
bration  data,  as  detailed  in  Ref.  1.  R^w  is  based  on  assumed 
local  warpage,  as  described  in  Section  3.4.  The  covariance 
of  temporal  grid  instability  is  related  to  lo  variations 

'V  ^ 

(An,  Aa ,  AC)  in  the  propagation  parameters  by 


tg 


=  H 


tg 


0,  o 

(An)2 


0 

0 


'V  o 

(Aa)2 


0 

0 

^  2 
(AC) 


H 


tg 


(  B-4 ) 


where  the  matrix  which  relates  parameter  variations  to  TD 
variations  is 


n 


( B-5 ) 


and  the  matrix  elements  are  defined  in  Apoendix  A  (Eqs.  A-6  to  A-8). 


The  diagonal  entries  of  P  are  denoted  by  p..  and  p^, 
and  the  error  quantity  in  common  usage  by  the  U.S.  Ccast  Guard 
(Ref.  3)  is 

2d  rms  error  =  2  VPN  +  PE  (B-6) 
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APPENDIX  C 

TABLES  OF  TD  VARIATIONS  AND  POSITION  ERRORS 


Tables  C-l  to  C-4  present  mean  TD  shifts  (TD  for 
the  season  minus  TD  for  September  1977)  and  random  (la)  TD 
variations  at  all  waypoints  for  four  seasons.  Table  C-5 
shows  the  magnitude  and  bearing  of  the  mean  position  error 
for  winter,  at  all  waypoints  for  the  station  combination  which 
yields  minimum  error  with  the  calibrated  grid  model.  Table 
C-6  presents  the  component  of  2d  rms  position  error  due  to 
random  seasonal  grid  instability,  at  all  waypoints  for  all 
three  station  combinations,  for  summer  and  winter.  Table  C-7 
presents  total  2d  rms  position  errors  at  all  waypoints  for 
all  three  station  combinations,  for  the  September  1977  data 
collection  period  and  for  winter. 
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TABLE  C-l 

PREDICTED  TD  VARIATIONS  AT  WAYPOINTS  FOR  SPRING 


T-2559 


WAYPOINT 

MEAN  TD  SHIFT 
( usee ) 

RANDOM  (la)  TD  VARIATION 
(usee ) 

TDX 

TDY 

TDZ 

TDX 

TDY 

TDZ 

A 

-0.001 

0.000 

-0.001 

0.  028 

0.015 

0.085 

B 

-0.001 

0.000 

-0.001 

0.032 

0.015 

0.044 

c 

-0.001 

0.000 

-0.001 

0.  024 

0.011 

0.041 

D 

0.000 

0.000 

ro.ooi 

0.023 

0  .011 

0.041 

E 

C.000 

0.000 

-0.  001 

0.020 

0.009 

0.04  3 

F 

0.  000 

0.  000 

-0*001 

0.013 

0.000 

0.043 

H 

0.000 

0. 000 

-0.001 

0.007 

0.003 

0.036 

i 

0.000 

0.000 

0.000 

0.001 

0.001 

0.027 

J 

0.000 

0.000 

0.000 

0.001 

0.000 

0.024 

K 

0.000 

0.000 

0.000 

0.002 

0*001 

0.016 

L 

0.000 

0.000 

0.000 

0.001 

0.000 

0.000 

M 

0.000 

0.000 

0.000 

0.001 

0.000 

0.009 

N 

0.000 

0.000 

0.000 

0.001 

0.000 

0.011 

0 

0.000 

0.000 

0.000 

0. 002 

0.002 

0.021 

p 

0.000 

0.000 

0.001 

0.007 

0.009 

0.041 

Q 

0.000 

0.000 

0.001 

0.002 

0.008 

0.056 

R 

0.000 

-0.001 

0.001 

0.010 

0.020 

0.057 

s 

0.000 

-0.001 

0.001 

0.010 

0.022 

0.060 

T 

-0.000 

0.000 

0.001 

0.007 

0.007 

0.031 

U 

-0.001 

-0.001 

0.001 

0.025 

0.019 

0.037 

V 

-0.001 

-0.001 

0.001 

0.028 

0.024 

o 

• 

o 

p 

<11 

w 

-0.001 

-0.001 

0.001 

0.037 

0.032 

0.048 

X 

-0.001 

-0.001 

0.001 

0.031 

0.038 

Y 

-0.001 

-0.001 

0.001 

0.033 

0.045 

0.061 

z 

-0.001 

-0.001 

0.001 

0.029 

0.059 

0.068 

AA 

-0.001 

-0.001 

0.001 

0.  030 

0.084 

0.068 

Bfl 

-0.001 

-0.002 

0.001 

0.042 

0.076 

0.064 

CC 

-0.001 

-0.002 

0.001 

0.044 

0.078 

0.064 

00 

-0.001 

-0.002 

0.001 

0.042 

0.080 

0.065 

EE 

-0.001 

-0.002 

0.002 

0.  034 

0.094 

0.069 

PF 

-0.001 

-0.002 

0.002 

0.  029 

0.1  08 

0.072 

ss 

-0.001 

-0.102 

0.002 

0.034 

0.110 

0.070 

HH 

•0.001 

-0. 103 

0.002 

0.030 

0.114 

0.069 

:  t 

-0.001 

-0.003 

0.002 

0.042 

0.123 

0.066 

THE  ANALYTIC  SCIENCES  CORPORATION 


TABLE  C-2 

PREDICTED  TD  VARIATIONS  AT  WAYPOINTS  FOR  SUMMER 

T-2560 


WAYPOINT 


MEAN  TD  SHIFT 

(.11  sec ) 

RANDOM 

lo)  TD  VARIATION 
(ysec ) 

TDX 

TDY 

TDZ 

TDX 

TDY 

TDZ 

0.000 

0.000 

0.000 

0.028 

Mi 

^4 

o 

• 

a 

0.065 

0.000 

0.000 

0.000 

0.C32 

0.015 

0.044 

o.coo 

0.000 

0.000 

0.024 

0.01  1 

0.041 

0.000 

0.000 

0.000 

o 

• 

o 

M 

III 

0.01  1 

0.041 

0.000 

0.  000 

0.  000 

0.020 

0.001 

0.C42 

0.000 

0.000 

0.000 

o 

• 

o 

0*005 

0.042 

0.000 

0.000 

0.000 

0.007 

0.00  ? 

0.036 

0.000 

0.000 

0.000 

0.001 

0.001 

0.027 

0.000 

0.000 

0.000 

0.001 

0.003 

0.024 

0.000 

0.000 

0.000 

0.0C2 

O.OC  1 

0.016 

0.000 

0.000 

0.000 

0.001 

0.003 

0.000 

0.000 

0.000 

0.000 

0.0C1 

0.003 

0.009 

0.000 

0.000 

0.000 

0.001 

0.003 

0.011 

0.000 

0.000 

0.000 

0.002 

0.003 

0.021 

0.000 

0.000 

0.000 

0.007 

0.003 

0.041 

0.000 

0.000 

0.000 

0.002 

0.003 

0.056 

0.000 

0.000 

0.000 

0.010 

0.02  ) 

0.057 

0.000 

0.000 

0.000 

0.010 

0.02  1 

0.060 

0.  000 

0.000 

0.000 

0.  000 

0.000 

0.000 

0.0  08 

0.025 

0.00  ’ 

0.01') 

■BjijB 

0.000 

0.000 

0.000 

0.028 

0.02<v 

0.045 

0.  000 

0.000 

0.000 

0.037 

0.03:: 

0.048 

0.000 

0.000 

0.000 

0.031 

0  •  03( 

0.056 

0.000 

0.  000 

0.000 

0.033 

0 .04? 

0.061 

0.000 

-0.001 

0.001 

0.029 

0.05‘> 

0.068 

0.000 

-0.001 

0.001 

0.  030 

0.054 

0.068 

0.000 

-0.001 

0.001 

0.042 

0.076 

0.064 

0.000 

-0.001 

0.  001 

0.  044 

0.07  5 

0.064 

0.000 

-0.001 

0.001 

0.042 

0.08  ) 

0.065 

0.000 

-0. 001 

0.001 

0.034 

o 

• 

o 

*> 

0.069 

0.000 

-0.001 

0.001 

0.029 

0.1  0( 

0.072 

0.000 

-0.00 1 

0.001 

0.034 

0  •  1  1  c 

0.C7C 

0.000 

-0.001 

0.001 

0.038 

0.11* 

0.069 

0.  000 

-0.001 

0.001 

0.042 

0.12' 

0.066 
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TABLE  C-3 

PREDICTED  TD  VARIATIONS  AT  WAYPOINTS  FOR  FALL 


T-2561 


WAYPOINT 

MEAN  TD  SHIFT 
(usee) 

RANDOM 

(la)  TD  VARIATION 
(usee) 

TDX 

TDY 

TDZ 

TDX 

TDY 

TDZ 

A 

0.000 

0.000 

0  .000 

0.020 

0.01S 

0.065 

B 

0.000 

0.000 

0.000 

0.032 

0.01S 

0.044 

C 

0.000 

0.000 

0.000 

0.024 

0.011 

0.041 

0 

0.000 

0.000 

0.000 

0.023 

0.011 

0.041 

e 

0.000 

0.000 

0.000 

0.020 

0.009 

0.043 

F 

0.000 

0.000 

0.000 

0.013 

0.006 

0.043 

H 

0.  000 

0.000 

0.000 

O.OC7 

0.003 

0.036 

i 

0.000 

0.000 

0.000 

0.001 

C.001 

0.027 

j 

0.000 

0.000 

0.000 

0.001 

0.000 

0.024 

K 

0.000 

0.000 

0.000 

0.002 

0.001 

0.016 

L 

0.000 

0.000 

0.000 

0.001 

0.000 

0.000 

M 

0.000 

0.000 

0.000 

0.001 

0.000 

0.009 

N 

0.000 

0.000 

0.000 

0.001 

0.000 

0.011 

o 

0.000 

0.000 

0.000 

0.002 

0.002 

0.021 

P 

0.000 

0.000 

0.000 

0.0C7 

0.009 

0.041 

0 

0.000 

0.000 

0.000 

0.002 

0.009 

R 

0.000 

0.000 

0.000 

0.010 

0.020 

0.057 

s 

0.000 

0.000 

0.000 

0.010 

0.022 

0.060 

T 

0.000 

0.000 

0.000 

O.OOB 

0.007 

0.032 

u 

0.000 

0.000 

0.000 

0.025 

0.019 

V 

0.000 

0.000 

0.000 

0 . 028 

0.024 

0.045 

V 

0.000 

0.000 

0.000 

0.037 

0.033 

0.048 

X 

0.000 

0.000 

0.00  0 

0.031 

0.036 

0.056 

Y 

0.000 

0.000 

0.000 

0.033 

0.045 

0.061 

z 

0.000 

0.  000 

0.000 

0.029 

0.059 

0.068 

AA 

0.000 

0.000 

0.000 

0.030 

0.064 

0.068 

BB 

0.000 

0*  000 

0.000 

0.042 

0.076 

o.oes 

CC 

0.000 

0.000 

0.000 

0.044 

0.078 

C.064 

OD 

0.000 

0.000 

0.000 

0.042 

o.oeo 

0.065 

EE 

0.000 

-0. 001 

0.000 

0.034 

0.094 

0.069 

FF 

0.000 

-0.001 

0.000 

0.029 

0.106 

0.072 

as 

0.000 

-0.001 

0.000 

0.034 

0.110 

0.070 

0.000 

-0.001 

0.000 

0.036 

0.1  IE 

0.069 

0.000 

-0.001 

0.000 

0.042 

0.124 

0.066 
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TABLE  C-4 

PREDICTED  TD  VARIATIONS  AT  WAYPOINTS  FOR  WINTER 


T-2562 


WAYPOINT 

MEAN  TD  SHIFT 
(usee) 

RANDOM  (lo)  TD  VARIATION 
(u«ec ) 

TDX 

TDY 

TDZ 

TDX 

TDY 

TDZ 

A 

-0.  027 

-0.117 

0.0*8 

0.026 

0.1 11 

B 

-0.027 

-0.079 

0.05* 

0.026 

0.076 

C 

-0.0*4 

-0.020 

-0.07* 

0.0*2 

0.019 

0.071 

D 

-0.0*1 

-0.019 

-0.07* 

0.039 

0.018 

0.071 

E 

-0.036 

-0.017 

-0. 078 

0.03* 

0.016 

0.07* 

F 

-0.023 

-0.010 

-0.077 

0.021 

0.010 

0.07* 

H 

-0.012 

-0.005 

-0.065 

0.011 

0.005 

0.062 

i 

-0.002 

-0*001 

-0.0*8 

0.002 

0.001 

0.0*6 

J 

-0.001 

0.000 

-0.0*2 

0.001 

0.000 

0.0*0 

K 

0.003 

0.002 

-0.029 

0.003 

0.002 

0.028 

L 

0.001 

0.001 

0.001 

O.OC1 

0.001 

0.001 

M 

0.002 

0.000 

0.017 

0.OC2 

0.000 

0.016 

N 

0.002 

0.000 

0.020 

0.002 

0.000 

0.019 

e 

-0 . 00* 

-0.00* 

0.038 

0.00* 

0.004 

0.036 

P 

-0.012 

-0.016 

0.073 

0.011 

0.015 

0.069 

a 

0.00* 

-0.015 

0.100 

0.00* 

0*014 

0.095 

R 

-0.019 

-0.035 

0.103 

0.018 

0.03* 

0.098 

s 

-0.019 

-0.0*0 

0.108 

0.  01B 

0.038 

0.103 

T 

-0.013 

-0.013 

0.056 

0.013 

0.012 

0.05* 

U 

-0*045 

-0.03* 

0. 066 

0.0*3 

0.032 

0.063 

V 

-0.050 

-0.0** 

0.080 

0.0*8 

0.0*2 

0.076 

» 

-0.066 

-0.058 

0.085 

0.063 

0.055 

0.081 

X 

-0.055 

-0.065 

0.100 

0.  053 

0.062 

111 

0> 

o 

• 

o 

Y 

-0. 059 

-0.051 

0.10B 

0.056 

0.077 

0.103 

2 

-0.052 

-0.106 

0.121 

0.0*9 

0.101 

0.115 

AA 

-0.055 

-0.11* 

0.121 

o 

• 

o 

M 

0.108 

0.116 

BB 

-o.ors 

-0.13* 

0.071 

0.129 

0.110 

CC 

-0.079 

-0. 1*0 

1 

0.075 

0.133 

0.1  09 

DD 

-0.076 

-0.1*3 

0.116 

0.072 

0.136 

0.111 

EE 

-0.062 

-0.168 

0.12* 

0.  059 

0.160 

0.1  18 

-0.052 

-0.191 

0.129 

O.CSO 

0.1  32 

0.123 

-0.061 

-0.197 

0.125 

0.  058 

0.1  87 

0.119 

-0.06* 

-0.205 

0.12* 

0.  061 

0.195 

0.118 

-0.07* 

-0.221 

0.119 

0.071 

0.211 

0.113 

C-5 


THE  ANALYTIC  SCIENCES  CORPORATION 


TABLE  C-5 

PREDICTED  MEAN  POSITION  ERRORS  AT  WAYPOINTS  FOR  WINTER 
DUE  ONLY  TO  VARIATIONS  IN  PROPAGATION  PARAMETERS 


T-2563 


WAYPOINT 

BEST 

TRIAD 

MEAN  WINTER  POSITION  ERROR 

MAGNITUDE  (ft) 

BEARING  FROM  NORTH 

A 

mxz 

200 

-66 

B 

MXZ 

92 

-85 

C 

UXZ 

64 

-82 

D 

uxz 

60 

-81 

E 

uxz 

58 

-77 

F 

MXZ 

45 

-66 

H 

UXZ 

34 

-55 

I 

uxz 

26 

-40 

J 

MXZ 

23 

-38 

K 

UXZ 

18 

-26 

L 

UXZ 

1 

45 

M 

MXZ 

9 

137 

N 

UXZ 

12 

138 

0 

uxz 

22 

125 

P 

UXY 

18 

159 

Q 

MXY 

22 

133 

H 

UXY 

30 

153 

S 

UXY 

33 

151 

T 

MXY 

14 

170 

U 

MXY 

31 

-176 

V 

UXY 

36 

178 

W 

UXY 

46 

179 

X 

MXY 

47 

171 

Y 

MXY 

56 

168 

Z 

MXY 

67 

162 

AA 

MXY 

72 

163 

SB 

MXY 

92 

171 

CC 

UXY 

97 

172 

DD 

UXY 

98 

171 

EE 

MXY 

107 

163 

FF 

MXY 

116 

158 

GG 

MXY 

129 

162 

HH 

UXY 

142 

163 

Cl 

MXY 

188 

166 

■‘Mean  position  errors  for  spring,  summer  and  fail  are 
less  than  1  ft. 
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TABLE  C-6 

PREDICTED  COMPONENT  OF  2d  RMS  POSITION  ERROR  AT 
WAYPOINTS  DUE  TO  RANDOM  GRID  INSTABILITY, 

FOR  SUMMER  AND  WINTER* 


2d  RMS  POSITION  ERROR  (It) 


WAYPOINT 

MXY 

MYZ 

MXZ 

SUMMER 

WINTER 

SUMMER 

WINTER 

SUMMER 

WINTER 

A 

67 

112 

301 

314 

223 

B 

j  57 

96 

126 

213 

103 

c 

41 

70 

86 

147 

71 

0 

36 

65 

82 

139 

6T 

E 

34 

56 

7B 

135 

64 

110 

F 

22 

37 

36 

1  00 

31 

86 

H 

12 

19 

41 

TO 

36 

65 

i 

3 

4 

27 

47 

29 

49 

J 

2 

2 

24 

41 

26 

44 

K 

4 

6 

17 

29 

20 

34 

«- 

1 

2 

1 

2 

1 

z 

H 

2 

4 

10 

17 

10 

18 

H 

3 

s 

12 

20 

13 

22 

0 

e 

13 

24 

41 

28 

48 

p 

20 

33 

31 

87 

62 

106 

0 

23 

42 

72 

122 

87 

148 

p 

34 

37 

BS 

144 

113 

192 

s 

37 

62 

92 

136 

126 

215 

T 

16 

27 

39 

66 

47 

81 

u 

34 

59 

36 

96 

74 

V 

41 

69 

71 

121 

98 

167 

w 

32 

88 

63 

142 

122 

208 

X 

53 

90 

100 

170 

132 

259 

r 

62 

106 

118 

201 

197 

336 

z 

73 

126 

131 

238 

266 

4  87 

AA 

60 

136 

157 

266 

309 

328 

BB 

103 

164 

279 

379 

647 

CC 

108 

iv  r 

282 

397 

677 

00 

109 

■5K. 

170 

290 

412 

702 

EE 

120 

204 

204 

346 

533 

909 

FF 

130 

222 

234 

399 

643 

1097 

6G 

144 

246 

229 

391 

679 

1157 

HH 

1 88 

270 

233 

397 

734 

1252 

it 

21  0 

356 

239 

407 

895 

1  525 

*Position  errors  for  spring  and  fall  are  within 
1  ft  of  those  for  summer. 
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TABLE  C-7 

PREDICTED  TOTAL  2d  RMS  POSITION  ERRORS  AT 
WAYPOINTS,  FOR  SEPTEMBER  1977  DATA 
COLLECTION  PERIOD  AND  WINTER 


2d  RMS  POSITION  ERROR  (ft) 


U7Z 


WINTER 


WINTER 

SEPT  1977 

855 

998 

540 

434 

352 

299 

321 

281 

307 

269 

253 

212 

195 

157 

SEPT  1977 

VINTER 

699 

764 

329 

358 

230 

250 

218 

237 

210 

228 

172 

185 

133 

143 

101 

108 

92 

99 

75 

80 

38 

38 

49 

51 

51 

54 

91 

99 

194 

212 

272 

297 

350 

383 

393 

430 

150 

164 

231 

253 

304 

333 

379 

415 

4  70 

515 

607 

665 

882 

966 

958 

1049 

1170 

1282 

1225 

1342 

1274 

1395 

1658 

1814 

2227 

2398 

2115 

2313 

2297 

2511 

2834 

3091 
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national  Radio  Consultative  Committee  (Geneva,  Switzerland), 


1974,  Vol .  V,  pp.  39-65. 
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